Antiphase magnetic boundaries in iron-based superconductors: A 
first-principle density-functional theory study 



Z. P. Yin and W. E. Pickett 
Department of Physics, University of California Davis, Davis, CA 95616 
(Dated: October 30, 2009) 

Superconductivity arises in the layered iron-pnictide compounds when magnetic long range order 
disappears. We use first principles density functional methods to study magnetic arrangements 
that may compete with long range order near the phase boundary. Specifically, we study the 
energetics and charge density distribution (through calculation of the electric field gradients) for 
ordered supercells with varying densities of antiphase magnetic boundaries. We quantify the amount 
by which Fe atoms with low spin moments at the antiphase boundaries have higher energies than Fe 
atoms with high spin moments away from the antiphase boundaries. These disruptions in magnetic 
order should be useful in accounting for experimental data such as electric field gradients and 
hyperfine fields on both Fe and As atoms. 



I. BACKGROUND AND MOTIVATION 

The discoveries of high temperature su- 
perconductivity in the FeAs-based layered 
compound^^^iViV . 11 . 12 . 13 of doped or 
compressed ZrCuSiAs-type i?FeAsO and XFeAsF, 
ThCr 2 Si 2 -type A/Fe 2 As 2 , and Cu 2 Sb-type AFeAs 
(R: rare-earth metal; X—Ca, Sr; M=Ca, Sr, Ba, 
Eu; A— hi, Na), with critical temperatures up 
to 56 K, has been attracting excitement in the 
condensed matter community. In the vicinity of 
room temperature, these compounds crystallize in 
tetragonal symmetry with no magnetic order. At 
some lower temperature (which can be in the range 
of 100 - 210 K), they undergo a first or second order 
phase transition to an orthorhombic structure and 
become antiferromagnetically ordered i 14 ' 15 ' 16 ' 17 The 
structural transition and magnetic order transition 
can happen simultaneously or successively depend- 
ing on the compound ] 14 ' 16 ! 17 It was confirmed both 
experimentally and theoretically that the magnetic 
order of Fe at low temperature is stripe-like antifer- 
romagnetism often referred to as spin density wave 
(SDW)i 14 ' 16 i 17 i 18 i 19 Upon doping or compressing, 
the magnetic order goes away and the materials 
become superconducting. 

Much theoretical work has been reported since 
the first discovery^ of LaFeAsOi_ x F x , with 
many aspects of these compounds having been 
addressed ^^^i^i^^i^i^ but with many ques- 
tions unresolved. A central question is what occurs 
at the SDW-to-SC (superconducting) phase transi- 
tion, and what drives this change, and more funda- 
mentally what microscopic pictures is most useful in 
this enterprise. In the i?FeAsO compounds, doping 
with carriers of either sign leads to this transition, 



even though there seems little that is special about 
the band-filling in the stoichiometric compounds. In 
the MFe 2 As 2 system, the SDW-to-SC transition can 
be driven with pressure (relatively modest, by re- 
search standards) without any doping whatever, ap- 
parently confirming that doping level is not an essen- 
tial control parameter. Some delicate characteristic 
seems to be involved, and one way of addressing the 
loss of magnetic order is to consider alternative types 
of magnetic order, and their energies. 

Many results, experimental and theoretical, indi- 
cate itinerant magnetism in this system, and LSDA 
calculations without strong interaction effects in- 
cluded correctly predict the type of antiferromag- 
netism observed. There is however the general 
feature that the calculated ordered moment of Fe 
is larger than the observed value. For example, 
neutron scattering experiment^ obtained the or- 
dered Fe magnetic moment of 0.36 [Ib in LaFeAsO, 
while calculation ;] 18 ! 21 ' 24 result in the much larger 
values 1.8-2.1 \xb- Neutron diffraction and neu- 
tron scattering experiment o 26 ' 27 ' 28 estimated the Fe 
magnetic moment in the SDW state of MFe 2 As 2 
(M=Ba, Sr, Ca) to be in the range of 0.8 /is to 
1.0 hb but our calculations (this work) give 1.6- 
1.9 jiB • This kind of (large) discrepancy of the 
ordered magnetic moment is unusual in Fe-based 
magnets, and there are efforts underway to under- 
stand the discrepancy as well as the mechanism un- 
derlying magnetic intereactions. 29 In addition, 57 Fe 
Mossbauer experiments 1 ^ ^^ 2 ^^ 4 ^!^ 6 . and 75 As 
NMR measurement a 16 ' 17 ' 37 further confirm the dis- 
agreement in magnetic moments and electric field 
gradients between experiments and ab initio calcu- 
lations in the SDW state. 

To explain these significant disagreements, it is 
likely that spin fluctuations in some guise play a 
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1,2,3: different As sites; A, B: different Fe sites. 

FIG. 1: The structure of FeAs layer in the Q-SDW state 
showing the antiphase boundary in the center of the fig- 
ure. Fe spin 1 (filled circle) and spin 2 (empty circle) 
have two different sites A ('bulklike') and B ('boundary- 
like'). As above Fe plane (filled square) and below Fe 
plane (empty square) have three sites 1, 2, 3 whose local 
environments differ. 



role in these compounds. The SDW instability is 
a common interpretation of the magnetic order in 
these compounds^ ' 14 ' 38 ' 39 which implicates the in- 
fluence of spin fluctuation in the magnetically dis- 
ordered state. An inelastic neutron scattering study 
on a single-crystal sample of BaFe2As2 by Matan et 
al^ showed anisotropic scattering around the an- 
tiferromagnetic wave vectors, suggestive of two di- 
mensional spin fluctuation in BaFe2As2- Such pos- 
sibilities must be reconciled with the existence of 
high energy spin excitations in the SDW state of 
BaFe2As2 as observed by Ewings et al£^. 

One of the simplest spin excitations is that aris- 
ing from antiphase boundaries in the SDW phase. 
Mazin and Johannes have introduced such "antipha- 
sons" and their dynamic fluctuations as being cen- 
tral for understanding the various phenomena ob- 
served in this class of materials^ The structural 
transition followed by the antiferromagnetic tran- 
sition, the change of slope and a peak in the dif- 
ferential resistivity dp(T)/dT at the phase transi- 
tions, and the invariance of the resistivity anisotropy 
over the entire temperature range can be qualita- 
tively understood in their scenario by considering 
dynamic antiphase boundaries (twinning of mag- 
netic domains) 42 

In this paper, we consider a class of magnetic ar- 
rangements derived from the stripe-like AFM phase: 
static periodic magnetic arrangements (SDWs) rep- 
resenting antiphase boundaries that require doubled, 
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quadrupled, and octupled supercells. We denote 
these orders as D-SDW, Q-SDW and O-SDW, re- 
spectively. Figure [1] shows the magnetic arrange- 
ments of Fe in the Q-SDW phase. Its unit cell 
is a 4x1x1 supercell of the SDW unit cell. An- 
tiphase boundaries occur at the edge and the center 
of its unit cell along a-axis (antiparallel/alternating 
Fe spins), the same as in the D-SDW and O-SDW 
states, the unit cells of which are 2x1x1 and 8x1x1 
supercells of the SDW unit cell, respectively. The 
D-SDW phase can also be viewed as a double stripe 
SDW phase. Based on the results of these states, 
we consider the effect of antiphase boundary spin 
fluctuations in explaining various experimental re- 
sults, which was discussed to some extent by Mazin 
and Johannes 42 The antiphase magnetic boundaries 
we consider here are the simplest possible and yet 
explain semiquantitatively many experimental re- 
sults by assuming that the dynamic average over 
antiphase magnetic boundaries can be modeled by 
averaging over several model antiphase boundaries. 
Our picture presumes that the antiphase boundary 
within the magnetically ordered state is in some 
sense a representative spin excitation in the FeAs- 
based compounds. 

The calculations are done using the linear aug- 
mented plane wave (LAPW) method as imple- 
mented in the Wien2k code^ 3 -, with both PW91 44 
and PBE4£ exchange-correlaton (XC) potentials. 
Several results have been double checked using the 
full potential local orbital (FPLO) i code4^41 with 
PW91 XC potential. 



II. LDA VS. GGA 

Whereas the local density approximation (LDA) 
for the exchange-correlation (XC) potential usu- 
ally obtains internal coordinates accurately, it has 
been foun d 21 ' 24 that LDA makes unusually large er- 
rors when predicting the As height z(As) in these 
compounds in either nonmagnetic (NM) or SDW 
states. The generalized gradient approximation 
(GGA) makes similar errors in the NM state, how- 
ever, GGA predicts very good values of z(As) in the 
SDW phase, as shown in Fig. [51 One drawback 
of GGA is that it enhances magnetis m 21 ' 24 in these 
compounds over the LDA prediction, which is al- 
ready too large compared to its observed value. For 
example, using experimental structural parameters, 
GGA (PBE) gives a Fe spin magnetic moment larger 
than LDA (PW91) by 0.3 \i B in the SDW state, and 
more than 0.6 fiB in the D-SDW state. The mag- 
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FIG. 2: The calculated errors of z(As) compared to ex- 
perimental values in the NM and SDW states when us- 
ing LDA (PW91) and GGA (PBE) XC functionals in 
CaFe 2 As 2 , SrFe 2 As 2 , BaFe 2 As 2 , and LaFeAsO. 



netic moment changes the charge density, roughly in 
proportion to the moment. For this reason, we have 
used Wien2K with PW91 (LDA) XC functional with 
its more reasonable moments to calculate the EFG 
and hyperfine field using experimental structural pa- 
rameters. We note that PBE and PW91 produce 
about the same EFG in the NM state. 



III. THE MAGNETIC MOMENT AND 
HYPERFINE FIELD OF IRON 

In the various antiphase boundary SDW states, 
the Fe atoms can assume two different characters: 
high-spin A site away from the antiphase boundaries 
and low-spin B site at the antiphase boundaries, as 
shown in Fig. [TJ The Fe atoms with the same site 
(A or B) in these states have about the same mag- 
netic moment and hyperfine field. For example, in 
BaFe2As2 in the static Q-SDW state, the spin mag- 
netic moment and hyperfine field for Fe at A site are 
1.59 \ib and 12.6 T, and for Fe at B site are 0.83 [is 
and 6.2 T. In the O-SDW state, the spin magnetic 
moments and hyperfine fields for the three (slightly 
different) A sites are 1.59, 1.60, 1.67 /is and 12.6, 
13.5, 1.37 T, respectively, and are 0.83 fiB and 6.1 
T for Fe at B site. 

As mentioned above, there are significant differ- 
ences of the ordered magnetic moment of Fe in 
the SDW state between calculated values and val- 
ues observed in neutron scattering (diffraction) ex- 
periments and/or Mossbauer experiments. Table U 



TABLE I: The experimental magnetic moment of Fe 
niFe (in unit of /xs) and the hyperfine field B^/ (in 
unit of Tesla) for Fe, and values calculated in the 
SDW and D-SDW ordered phases, using Wien2K with 
PW91 for the MFe 2 As 2 (M=Ba, Sr, Ca), LaFeAsO 
and SrFeAsF compounds. The experimental values are 
in all cases much closer to the D-SDW values (with 
its maximally dense antiphase boundaries) than to the 
SDW values^^^i^.^^^ For the Fe magnetic mo- 
ment, results from both FPLO (denoted as FP) and 
Wien2k (denoted as WK) are given. Because these meth- 
ods (and other methods) differ somewhat in their assign- 
ment of the moment to an Fe atom, the difference gives 
some indication of how strictly a value should be pre- 
sumed. 
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compound 
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Bhf 




B/i/ 












FP 


WK 


WK 


FP 


WK 


WK 


BaFe 2 As 2 


0.8 


5.5 


1.78 


1.65 


13.6 


0.90 


0.80 


5.4 


SrFe 2 As 2 


0.94 


8.9 


1.80 


1.68 


13.9 


0.98 


0.91 


6.1 


CaFe 2 As 2 


0.8 




1.63 


1.53 


12.4 


0.77 


0.71 


4.7 


LaFeAsO 


0.36 




1.87 


1.77 


14.9 


0.50 


0.48 


3.6 


SrFeAsF 




4.8 




1.66 


14.5 




0.32 


2.3 



shows Fe spin magnetic moment and hyperfine field 
in the SDW and static D-SDW state calculated by 
FPL07 and Wien2K using PW91 XC functional. 



IV. ENERGY DIFFERENCES 

Since Fe atoms in all these antiphase boundary 
SDW states have basically two spin states (high-spin 
state at A sites and low-spin state at B sites), in a 
local moment picture one might expect that energy 
differences could be related to just the two corre- 
sponding energies. From our comparison of energies 
we have found that this picture gives a useful ac- 
count of the energetic differences. 

The SDW and the D-SDW state are the simple 
cells in this regard. The former doesn't have any 
low spin Fe and the latter doesn't have any high 
spin Fe, so these two define the high spin (low en- 
ergy) and low spin (high energy) "states" of the Fe 
atom. Table [II] shows the total energies per Fe (in 
meV) of the magnetic phases compared to the non- 
magnetic stater, for BaFe2As2, SrFe2As2, CaFe2As2, 
LaFeAsO and SrFeAsF. The high spin and low spin 
energies vary from system to system. The energy of 
the Q-SDW state has also been calculated, and it 
can be compared with the average of high spin and 
low spin moment energies (last column in Table [TTJ) . 
The reasonable agreement indicates that corrections 



TABLE II: Calculated total energies (meV/Fe) com- 
pared to NM state of the various SDW states (SDW, 
D-SDW, Q-SDW) in the MFe 2 As 2 (M=Ba, Sr, Ca), 
LaFeAsO and SrFeAsF compounds. The energy tabu- 
lated in the last column, labeled Q', is the average of the 
high spin (SDW) and low spin (D-SDW) energies, il- 
lustrating that the energy of the Q-SDW ordered phase 
follows this average reasonably well. The level of agree- 
ment indicates to what degree 'high spin' and 'low spin' 
is a reasonable picture of the energetics at an antiphase 
boundary. 



compound 


SDW 


D 


Q 


Q' 


BaFe 2 As 2 


-73 


-6 


-36 


-39 


SrFe 2 As 2 


-91 


-11 


-46 


-51 


CaFe 2 As 2 


-66 


-8 


-33 


-37 


LaFeAsO 


-143 


-61 


-94 


-102 


SrFeAsF 


-73 





-40 


-37 



beyond this simple picture are minor. 

The energy cost to create an antiphase boundary 
is (roughly) simply the cost of two extra low spin Fc 
atoms versus the high spin that would result without 
the antiphase boundary. This difference is found to 
vary by over a factor of two, in the range of 40- 
90 meV per Fe for this set of five compounds. The 
reason for the variation is not apparent; for example, 
it is not directly proportional to the Fe moment (or 
its square). 

Analogous calculations were also carried out for 
the large O-SDW cell for BaFe 2 As 2 . As for the 
other compounds and antiphase supercells, the Fe 
moments could be characterized by a low spin atom 
at the boundary and high spin Fe elsewhere. The 
energy could also be accounted for similarly, analo- 
gously to Table III 



V. THE ELECTRIC FIELD GRADIENT 

The EFG at the nucleus of an atom is sensitive 
to the anisotropy of the electron charge distribution 
around the atom. A magnitude and/or symmetry 
change of the EFG implies the local environment 
around the atom changes, which can be caused by 
changes in bonding, structure, or magnetic ordering. 
In BaFe2As2 during the simultaneously structure 
transition from tetragonal to orthorhombic and mag- 
netic order transition from non-magnetic to SDW 
order at about 135 K, the EFG component V c along 
the crystal c-axis drops rapidly by 10% and the 
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the largest component V zz is changed from along c- 
axis to in the ab plane— The abrupt EFG change 
reflects a large change in the electron charge distri- 
bution around As sites and highly anisotropic charge 
distribution in the ab plane. A similar thing hap- 
pens in CaFe2As2 except that the V c component 
of CaFe2As2 in the non-magnetic state is five times 
that in BaFe2As2, and doubles its value at the struc- 
tural and magnetic transition at 167 K when it 
goes to the SDW phased The different behaviour 
of the EFG change across the phase transition in 
BaFe2As2 and CaFe2As2 may be due to the out-of- 
plane alkaline-earth atom (Ba and Ca in this case), 
which influences the charge distribution around As 
atoms. It also indicates that three dimensionality 
is more important in A/Fe2As2 than in i?FeAsO, 
which is evident in the layer distance of the FeAs 
layers reflected in the c lattice constant of these 
compounds. The c lattice constant of CaFe 2 As2 is 
significantly smaller than BaFe2As2, so that the in- 
terlayer interaction of the FeAs layers is stronger 
in CaFe2As2, therefore the charge distribution in 
CaFe2As2 is more three dimensional like than in 
BaFe2As2, which can be clearly seen in their Fermi 
surfaces (not shown). 



A. EFG of As atoms 

The EFG of As can be obtained from the 
quadrupolc frequency in nuclear quadrupolar reso- 
nance (NQR) measurement in nuclear magnetic res- 
onance (NMR) experiment. The NQR frequency can 
be written as 



V Q 



ZeQV zz {l + if/?>) 1 ' 2 
21(21 -l)h 



(1) 



where Q (~ 0.3b, b=10~ 24 cm 2 ) is the 75 As 
quadrupolar moment, V zz is the zz component of 
As EFG, 
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(2) 



to larger than one, indicating the principle axis for 



is the asymmetry parameter of the EFG, 1=3/2 is 
the 75 As nuclear spin, h is the Planck constant. In 
LaO .gF .iFeAs, Grafe et alML reported ^g=10.9 
MHz, and r?=0.1, which gives V zz ~ 3.00 x 10 21 
V/m 2 . (Note: for the value of EFG, the unit 10 21 
V/m 2 is commonly used, and we will adopt this unit 
for all EFG values below) The experimental value 
3.0 agrees satisfactorily with our result 21 of 2.7 cal- 
culated by Wien2K code. 



In BaFe2As2 and CaFe2As2, NMR experiments 
suggest the V c component of the EFGs of As are 
0.83 and 3.39 respectively at high temperature in 
the NM stateS ) 16 i 17 while our calculated values are 
1.02 and 2.35, respectively. The difference is in the 
right direction and right order of magnitude though 
not quantitatively accurate. However, the EFGs cal- 
culated in the SDW state don't match experimental 
observations at all. In BaFe2As2 from 135 K down to 
very low temperature, V c remains around 0.62 and 
7] = ^ Va \yJ b ^ changes from 0.9 to 1.2. Our calculated 
results in the SDW state gives V Q =1.34, V fc =-1.47, 
and V c =-0.13, which gives r] ss 20. The calculated 
results substantially underestimate V c and overem- 
phasize the anisotropy in the ab plane. 

We now consider whether these discrepancies can 
be clarified if antiphase boundary is considered. In 
the static D-SDW, Q-SDW and O-SDW states, the 
surrounding environment of As sites change. De- 
pending on the magnitudes (high spin or low spin) 
and directions (parallel or antiparallcl) of the spins 
of their nearest and next nearest neighboring Fe 
atoms, As atoms generally have three different sites 
(1, 2, and 3) as shown in Fig. [TJ In the static D- 
SDW state, As atoms have similar site 1' and 3'. 
As shown in Table IIIH the calculated quantities for 
these states cannot directly explain the experimental 
observed values neither. 

However, they may be understandable if the an- 
tiphase boundary is dynamic, i.e., any As atom in 
a given measurement can change from site 1 to site 
2 and/or site 3, when its nearby Fe atoms flip their 
spin directions. These time fluctuations have be rep- 
resented by an average over configurations, and we 
consider briefly what arises from a configuration av- 
erage over our SDW phase and three short-period 
ordered cells (having different densities of anti-phase 
boundaries). In the O-SDW state, for example, if an 
As samples 25%, 25% and 50% "time" as As sites 1, 
2 and 3, respectively during the experimental mea- 
surement, then the expectation values are V a =0.99, 
V 6 =-0.29, V c =-0.71, 77=1.80, B^l.43 T. This sim- 
ple consideration already match much better with 
experimental^ observed V c ~ 0.62, Hj n =1.4 T, ex- 
cept r\ which is around 1.2. The actual situation 
could be much more complicated, being an aver- 
age over all the sites in all the static D-SDW, Q- 
SDW and O-SDW states, and other more compli- 
cated states. Considering the relatively small dif- 
ferences of the EFGs at the same site for Q-SDW 
and O-SDW order, As sites in other static antiphase 
boundary SDW states should be able to be classi- 
fied to site 1, 2 and 3 as in the Q-SDW and O-SDW 



TABLE III: The calculated EFG component V„, V b , V c 
(in unit of 10 21 V/m 2 ) , the asymmetry parameter r/, 
spin magnetic moment of As (hb), hyperfine field at the 
As nuclei (Tesla) of BaFe 2 As 2 in the SDW, D-SDW, Q- 
SDW and O-SDW states. Experimentally, V c is around 
0.62, r\ is in the range of 0.9 to 1.2, and the internal field 
at As site parallel to c axis is about 1.4 TM See text for 
notation. 



state 


site 


Va 


Vb 


V c 


V 


Bhf 


SDW 


1 


1.21 


-1.32 


0.11 


23.0 





D-SDW 


1' 


0.63 


0.07 


-0.70 


0.8 







3' 


0.66 


0.37 


-1.03 


0.28 


2.1 


Q-SDW 


1 


1.17 


-1.38 


0.21 


12.1 







2 


0.99 


-0.72 


-0.27 


6.33 


1.0 




3 


0.90 


0.51 


-1.41 


0.28 


2.1 


O-SDW 


1 


1.25 


-1.44 


0.19 


14.2 







2 


0.98 


-0.71 


-0.27 


6.26 


1.1 




3 


0.87 


0.50 


-1.37 


0.27 


2.3 



states. 



B. EFG of Fe atoms 

The EFG of Fe can be obtained from the elec- 
tric quadrupole splitting parameter derived from 
Mossbauer measurements. The electric quadrupole 
splitting parameter can be written as 

. _ 3eQV«(l + 7, 2 /3) 1/2 



which equals E 7 AEq/c, where Q ~ 0.16b is the 57 Fe 
quadrupolar moment, E 7 is the energy of the 7 ray 
emitted by the 57 Co/Rh source, AEq is the elec- 
tric quadrupole splitting parameter from Mossbauer 
data given in the unit of speed, and c is the speed 
of light in vacuum. By fitting the Mossbauer spec- 
tra, one also obtains the isomer shift S and average 
hyperfine field B^/. 

We consider whether the dynamic antiphase 
boundary spin fluctuation picture can also clarify the 
comparison between calculated and observed EFG 
of Fe. We take SrFe2As2 as an example. As shown 
in Table IIV1 in the NM state, the calculated value 
Vq=0.98 agrees well with the Vq w 0.83 at room 
temperature. Vq calculated in the SDW state (0.68) 
agrees rather well with experimental value about 
0.58 at 4.2 K. Vq calculated in the D-SDW state 
(0.80) is somewhat larger than that in the SDW 
state, but not . Regarding EFG, the biggest dif- 
ference between SDW and D-SDW is the asymme- 
try parameter-it is 0.61 in the SDW and only 0.11 



TABLE IV: The calculated EFG component V , V 6 , 
V c (in unit of 10 21 V/m 2 ) , the asymmetry parame- 
ter n=\V xx -'Vyy\/\'Vzz\ (here jV zz | > \V XX \ and |Vj,j,|), 
Vq = |V zz |/(1+77 2 /3) 1/2 of Fe in SrFe 2 As 2 in the NM, 
SDW, D-SDW, Q-SDW and O-SDW states. Experimen- 
tally, Vq is around 0.83 at room temperature in the non- 
magnetic state, and it is about 0.58 at 4.2 K.— 



state 


site 




v b 


v c 


V 


Vq 


NM 





-0.49 


-0.49 


0.98 





0.98 


SDW 


A 


-0.58 


-0.14 


0.72 


0.61 


0.68 


D-SDW 


B 


-0.30 


-0.51 


0.81 


0.26 


0.80 


Q-SDW 


A 


-0.63 


-0.06 


0.69 


0.83 


0.62 




B 


-0.56 


-0.49 


1.05 


0.07 


1.05 



in the latter. Further experiments are required to 
clarify this difference. 

VI. SUMMARY 

Experiments generally indicate that an itinerant 
magnetic moment, magnetic (SDW) instability, and 
spin fluctuations are common features of the Fe- 
based superconductors. In this paper, we have 
studied the energetics, charge density distribution 
(through calculation of the electric field gradients, 
hyperfine fields and magnetic moments) for ordered 
supercells with varying densities of antiphase mag- 
netic boundaries, namely the SDW, D-SDW, Q- 
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SDW and (for very limitied cases) O-SDW phases. 
Supposing dynamic magnetic antiphase boundaries 
are present, and that the spectroscopic experiments 
average over them, we can begin to clarify several 
seemingly contradictory experimental and computa- 
tional results. 

Our calculations tend to support the idea that 
antiphase boundary magnetic configurations can be 
important in understanding data. The fact that 
the decrease in moment is confined to the antiphase 
boundary Fe atom does not mean that a local mo- 
ment picture is appropriate; in fact, exactly this 
same type of local spin density calculations provide 
a description of magnetic interaction that is at odds 
with a local moment pictured The calculated en- 
ergy cost to create an antiphase boundary is however 
rather high for the cases we have considered, and 
this value would seem to restrict formation of an- 
tiphase boundaries at temperatures of interest. Cal- 
culations that treat actual disorder, and dynamics 
as well, would be very helpful in furthering under- 
standing in this area. 
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